Phase shift experiments identifying Kramers doublets in a chaotic superconducting 

microwave billiard of threefold symmetry 
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The spectral properties of a two-dimensional microwave billiard showing threefold symmetry have 
been studied with a new experimental technique. This method is based on the behavior of the 
eigenmodes under variation of a phase shift between two input channels, which strongly depends 
on the symmetries of the eigenfunctions. Thereby a complete set of 108 Kramers doublets has been 
identified by a simple and purely experimental method. This set clearly shows Gaussian unitary 
ensemble statistics, although the system is time-reversal invariant. 
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In a recent experiment with a superconducting mi- 
crowave resonator, which possesses threefold symmetry 
(C3 symmetry), first signatures for spectral properties 
according to Gaussian unitary ensemble (GUE) statis- 
tics have been observed for a time-reversal invariant sys- 
tem 0. Discrepancies between experiment and theory 
remained, although these results confirmed the theoret- 
ical prediction 2], that time-reversal invariant systems 
may show properties, characteristic of violation of time- 
reveral invariance (TRI). A comparison with numerical 
results allowed to understand the discrepancies, which 
were due to a less than perfect identification of Kramers 
doublets that are supposed to display GUE statistics. 
The only criterion used to identify these doublets, which 
are split due to C3 symmetry breaking perturbations, was 
the spacing between consecutive resonances. The results 
of Q were even worse and demanded for an explanation 
by ambitious numerical studies. In the present paper we 
report on a new experimental technique, which identi- 
fies the doublets properly and uniquely, and agreement 
between experiment and theory is obtained within sta- 
tistical uncertainty. 

Billiards are model systems, which allow to study the 
properties of classically chaotic dynamics: a point-like 
particle moving without losses inside a closed bound- 
ary, where it is reflected specularly at impact, possesses 
non-integrable equations of motion, if the boundary is 
suitably shaped (see, e.g., 0). These systems have a 
quantum mechanical analogue, where the chaotic behav- 
ior manifests itself in the properties of the eigenfunctions 
and spectra It is common knowledge that the statis- 
tical fluctuations of the eigenvalues of a quantum system 
with a classically chaotic counterpart can be described 
by random matrix theory 0, . For chaotic quantum 
systems one usually finds fluctuations similar to those 
of the Gaussian orthogonal ensemble (GOE), which re- 
flects TRI. It has been recognized, however, that certain 
systems with TRI show spectral fluctuations which fol- 
low the GUE and thus mimic a violation of TRI |2|. A 
simple example for a time-reversal invariant system with 



GUE statistics is a two-dimensional chaotic billiard, the 
highest symmetry of which is C3. Such a billiard has 
been used in pj and in the present work (see fig.0. As 
was pointed out in Q, Q the Hamiltonian of this triangu- 
lar billiard possesses three classes of eigenfunctions 
with different symmetry properties (I = — 1,0, +1). A 
rotation R of 120° transforms these eigenfunctions ac- 
cording to = exp(27rzi/3)<3> (i) . Two of the classes 
of wave functions are degenerate due to TRI and one 
finds a spectrum composed of singlet (/ = 0) and doublet 
(I = ±1) modes. The doublets were predicted to show 
GUE statistics, while the singlets should follow GOE 
statistics. In a complex representation time-reversal will 
transform one eigenfunction of a doublet into its com- 
plex conjugate and thus interchanges the two eigenfunc- 
tions, i.e. = T*( +1 ) = (*( +1 ))*. This is an 
example for Kramers degeneracy: whenever TRI holds, 
i.e. [H,T] — 0, the time-reversed quantity TvP of an 
eigenfunction is also an eigenfunction of H. Thus, if 
* ^ T*, TRI forces a degeneracy (see, e.g., Q). This 
is known as Kramers theorem and its consequences play 
an important role in studies of systems with half-integer 
spin [<j. 

Two-dimensional quantum billiards can be studied ex- 
perimentally with the help of macroscopic analog systems 
|l0l lll| : in flat cylindric resonators only transverse mag- 
netic modes exist below a critical frequency. These modes 
have electric field vectors perpendicular to the bottom 
and lid of the resonator, while the magnetic field lines lie 
parallel to them. The electric field E n (r) — <p n (x,y)e z 
as well as the magnetic field is fully described by the 
solutions 4> n (x,y) of the scalar Helmholtz equation 
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with Dirichlet boundary conditions. There is a complete 
analogy between the two-dimensional Helmholtz equa- 
tion for a flat microwave cavity and the Schrodinger equa- 
tion describing a quantum billiard of the same shape. 
The wave numbers k n of the electromagnetic resonator 
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correspond to the energy eigenvalues e n of the Hamilto- 
nian of the quantum billiard, i.e. fc^ oc e n . 

The usual way to determine the eigenfrequencies of a 
microwave cavity /„ = fc„c/(27r), where c denotes the 
speed of light, is to measure the power transmitted from 
one antenna emitting microwaves into the cavity to an- 
other antenna receiving the microwaves. As a result one 
gets a resonance line (see fig. EJ, where the positions of 
the maxima correspond to the eigenfrequencies /„, while 
the widths of the resonances Af n reflect the damping, 
characterized by the quality factor Q = f n /Af n . The 
widths of the resonances define the experimental limit of 
resolution. Thus, a superconducting cavity of high Q is 
a prerequisite for resolving small splittings of degenerate 
eigenvalues, which stem from mechanical imperfections, 
and for identifying Kramers doublets properly. In our ex- 
periment a resonator made from lead plated copper was 
used. The measurements were carried out in a cryostat 
with liquid helium at a temperature of 4.2 K. The res- 
onator becomes superconducting at temperatures below 
approximately 7.2 K and possesses a quality factor of the 
order of 10 4 . The microwave signals are processed by a 
HP-8510C network analyzer, which also controls the sig- 
nal generator providing the microwaves. This technique 
has been used in several experiments |lOl 1 1 1 1 | , including 
the first measurements with the cavity under considera- 
tion 0, and now was modified to improve the results in 
the following way (see fig-CJ. 
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FIG. 2: Part of a transmission spectrum measured at T = 
4.2 K. The maxima of the resonance line mark the resonance 
frequencies. One clearly observes a signal-to-noise ratio of 
about four orders of magnitude (40 dB). Singlets and doublets 
are marked by S and D, resp. 



while the output power is coupled out through a third 
channel. The positions of the three antennas were chosen 
in a way, which preserves the threefold symmetry of the 
billiard. A Midwest Microwave PWD-2533-02-SMA-79 
power divider splits the microwave power provided by 
the signal generator into two signals. These two input 
signals have a certain phase shift 
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FIG. 1: Schematic view of the experimental setup. The 
chaotic microwave cavity (left side) possesses threefold sym- 
metry. The microwaves provided by port 1 (PI) of the net- 
work analyzer are split into two signals before they are cou- 
pled into the cavity by two different channels. The phase 
shift A<E> between these two channels can be varied by a phase 
shifter. Power is coupled out by a third antenna and led back 
to the network analyzer, where the transmission to port 2 
(P2) is measured. 

An additional antenna is used for feeding microwave 
power into the resonator. Thus, power is coupled into the 
system through two different channels at the same time, 



where A$ can be varied by an ARRA 9428A phase 
shifter, while A<I>o is an offset, which is constant for a 
given frequency and is due to the different lengths of the 
two signal paths. The phase shifter and the power divider 
operate up to a frequency of 18 GHz. For frequencies 
below 6 GHz the maximum achievable A$ is less than 
360°. For the frequency range from 8 GHz up to 18 GHz 
50 different transmission spectra were measured, cover- 
ing phase shifts from A$ w 0° up to A<1> « 350° in steps 
of 7°. Data points were taken in steps of 100 kHz, i.e. a 
resolution of the order of 10 -5 , which has to be compared 
with 1/Q « 10~ 4 . Furthermore, conventional transmis- 
sion spectra - as described above - have been recorded 
in the frequency range from to 18 GHz with a spectral 
resolution of 50 kHz (see fig. EJ) ■ This was done for two 
different antenna combinations, with a fourth antenna 
not used for the phase shift experiment. These transmis- 
sion spectra served to identify the resonance frequencies, 
while the results of the phase shift measurements were 
utilized to identify the eigenmodes as singlets or dou- 
blets. 

As we had conjectured, singlets and doublets typically 
show a very different behavior, when the phase shifter 
setting A$ is varied. For a typical singlet mode (fig. |3J) 
the shape of the resonance line does not change, although 
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FIG. 3: Resonance line for a typical singlet. The shape of 
the line depends only weakly on the phase shift varied. Nev- 
ertheless, there is a change in amplitude of more than one 
order of magnitude, i.e. a difference of more than 10 dB in 
the logarithmic scale. 

there is a change in the amplitude: the transmitted power 
at resonance depends on A<3> and shows a minimum and 
a maximum, corresponding to destructive or constructive 
interferences. Because of the special symmetries of the 
antennas and of the singlet wave functions these should 
occur at a distance of 180°. However, the symmetries are 
not perfect and thus the minimum amplitude is non zero, 
but for some singlets the intensities at the maximum and 
the minimum differ by two orders of magnitude. Further- 
more minimum and maximum do not always have an ex- 
act phase difference of 180° as expected, but 179° ± 33°. 
There is a correlation between the deviation from 180° 
and the minimum intensity. Large shifts only appear, 
when the minimum and the maximum intensity differ 
only marginally. The shape of the resonance, however, 
does not change due to those effects. 

In contrast to the singlet levels, typical doublets show 
a strong variation of the shape of the resonance curve 
(fig-EI- The relative and the absolute height of each of 
the two peaks within a doublet depends on the phase 
difference and sometimes an almost complete vanishing 
of one of the peaks is achieved. In some cases a dou- 
blet can only be identified for certain values of the phase 
shift. However, the behavior of the doublet amplitudes 
depending on A$ is much harder to understand, as one 
has to deal with two eigenmodes being excited simul- 
taneously with different amplitudes, the ratio of which 
depends on the properties of the two individual states 
involved and is uncontrollable in our setup. Neverthe- 



FIG. 4: Resonance line for a typical doublet. The shape of 
the line strongly depends on the phase shift varied. Note that 
this is a semi- logarithmic plot, which compresses the changes 
in amplitude. 

less, a strong influence of the induced phase shift on the 
shape of a resonance line composed of a doublet is obvi- 
ous (see fig. 0J), which provides an important additional 
criterion for distinguishing doublets from singlets. 

With the help of the method described above, 102 sin- 
glets and 108 doublets have been identified unambigously 
in the frequency range up to 18 GHz. The sets of eigen- 
frequencies /, were unfolded using Weyl's law 

Xi := N Wey i(fi) = vijf + v 2 fi + v 3 (3) 

This procedure leads to spectra Xi with normalized mean 
level spacings of 0.988 and 1.006 for singlets and doublets, 
respectively. As the number of levels in the two spectra is 
only about 100, we chose to analyze the long-range corre- 
lations of the spectra. Evaluating the Dyson-Mehta (A3) 
statistics yields perfect GUE behavior for the dou- 
blets (see fig-EJ), including saturation 0,0], i.e. A 3 (L) 
stops increasing for L > L max . The value of L max w 10 
is consistent with an estimate of the length of the short- 
est peridic orbit in the billiard. All this indicates that 
the corresponding spectrum is indeed complete. Other- 
wise deviations would occur, as observed in our previous 
work For the singlets we still observe deviations from 
GOE for small spacings (see fig. |SJ . As there are fewer 
singlets than doublets, this has to be due to missing lev- 
els. A check with a newly developed method based upon 
the assumptions that resonances below a certain critical 
strength are not detected and that the positions of the 
missing levels are distributed randomly indeed indicates 
that we lost about 5 percent of the singlets, while we have 
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all of the doublets [15J . It is most likely that some of the 
singlet modes cannot be detected with our antennas, be- 
cause their positions show the same symmetry properties 
(C3) as the singlet wave functions. We cannot circumvent 
this problem without destroying the C3 symmetry, e.g. 
when inducing a perturbing body. Thus, if one antenna 
lies on a nodal line, all antennas do so at frequencies cor- 
responding to singlets, while this is not the case for the 
doublets. These considerations are also supported by re- 
sults we obtain for the Weyl coefficients v\, 1)2, and i>3, 
which were determined by a fit of eq. © to the exper- 
imental data. The coefficients v\ and 1*2 are connected 
to area A and perimeter U, resp., of the billiard. We 
found A = 0.03 m 2 and U = 1.16 m for the singlet 
spectrum and A — 0.03 m 2 and U = 0.81 m for the 
doublets. Therefore, the Weyl parameters of the doublet 
spectrum are much closer to the the physical dimensions 
of the cavity, which are approximately A = 0.03 m 2 
and U = 0.82 m at room temperature. Note, that 
for Dirichlet boundary conditions a larger perimeter U 
corresponds to a smaller number of eigenfrequencies ^3] • 

We can be sure, however, not to have spurious singlet 
levels. This was the case in jjj, where we did not have 
a reliable method of separating the subspectra, which 
caused misassignments, i.e. lost and extra eigenvalues 
in the resulting spectra of 196 singlets and 213 doublets 
up to 25 GHz. Thus, these spectra showed deviations 
from GOE- and GUE-like statistics, respectively 0. Also 
there was no saturation behavior of the Dyson-Mehta 
statistics (see fig- EJ - Now saturation is obtained both 
for singlets and doublets. Although the modified setup 
restricts us to a smaller frequency range, it nevertheless 
allows a much more reliable identification and classifica- 
tion of states as singlets or doublets. 

To conclude, the spectrum of a microwave cavity 
with C3 symmetry shows degenerate and non-degenerate 
modes. Under variation of a phase shift between two 
feeding antennas these classes of states behave differently. 
Therefore, the experimental technique presented allows 
us for the first time to classify singlets and doublets as 
such. For the doublet spectrum fluctuation properties of 
the eigenvalues perfectly matching GUE statistics have 
been observed for a system with TRI - with 108 Kramers 
doublets only. This constitutes a clear improvement upon 
earlier results, where additional numerical information 
was needed to fully identify doublets |l| or to explain the 
effects of level splitting and level loss |3J . 
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FIG. 5: Long range correlations for the singlets and doublets: 
in the upper half the Dyson-Mehta statistics obtained in the 
framework of our previous work is shown. As one can see 
there are distinct deviations from GOE for the singlets. Also 
there is no saturation behavior for singlets and doublets. In 
the lower part the Dyson-Mehta statistics for the results of 
the present work are plotted. The statistics show for both 
singlets and doublets saturation, although the singlets still 
deviate slightly from GOE for small interval lengths L. 
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